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A novel approach to solve the Faddeev equation for three-body scattering at arbitrary energies 

is proposed. This approach disentangles the complicated singularity structure of the free three- 

nucleon propagator leading to the moving and logarithmic singularities in standard treatments. 

The Faddeev equation is formulated in momentum space and directly solved in terms of momentum 

vectors without employing a partial wave decomposition. In its simplest form the Faddeev equation 

for identical bosons, which we are using , is an integral equation in five variables, magnitudes of 

00 ' relative momenta and angles. The singularities of the free propagator and the deuteron propagator 

C^ , are now both simple poles in two different momentum variables, and thus can both be integrated 

^^ . with standard techniques. 

(N 

0\ \ I. INTRODUCTION 



O 



PACS numbers: 21.45.-v,25.10.+s 



r^ , In 1960 L.D. Faddeev formulated his mathematically rigorous scattering theory for three particles by proposing a 

set of three coupled integral equations, which do have a unique solution [l|. In the first numerical realizations of 
this approach, separable interactions were introduced, reducing the three-body equations to a set of one-dimensional 
^ ' coupled integral equations, whose numerical solutions was feasible at the time [3, Q- Despite a simplification due 
P^ \ to the choice of the two-body interaction, the standard formulation of the momentum space Faddeev equations in 
the continuum contains for the free three-body propagator a complicated singularity structure within the integral 
kernel 0, H, 0] • These complications arise because the position of the propagator cut does not only depend on the 
total energy of the system but also on external momentum variables on a grid (thus the term moving singularities). 
In addition branch points occur leading to logarithmic singularities. 

^s) , It is possible to carry out the integration along a path in the complex plane, thus avoiding to directly deal with these 

^^ ' singularities, but indirectly imposing conditions on the analytical properties of the two-body force. This method of 
contour deformation was introduced for separable potentials by Hetherington and Schick Q and perfected by Cahill 

r — . and Sloan (8[. The need to use realistic forces, which are predominantly local, led to methods integrating the 

C^ ' singularities on the real momentum axis @- It took until the 1980's until the Faddeev equations were solved in the 

00 . continuum with a realistic nucleon-nucleon (NN) force as input [l^, llH ■ 

^^ ' During the last two decades calculations of nucleon-deuteron scattering experienced large improvements and refine- 

^ , mcnts. It is fair to say that below about the pion production projectile energy the momentum space Faddeev equations 

k>( ' ^*-"' three-nucleon scattering can now be solved with high accuracy for the most modern two- and three-nucleon forces. 

rN ' A summary of these achievements can be found in Refs. [l3, llJ, [iJ, Ej 13 ■ The approach described there is based on 
S ■ using angular momentum eigenstates for the two- and three-body systems. This partial wave decomposition replaces 
■ " " ' the continuous angle variables by discrete orbital angular momentum quantum numbers, and thus reduces the number 
of continuous variables to be discretized in a numerical treatment to two. For low projectile energies the procedure 
of considering orbital angular momentum components appears physically justified due to arguments related to the 
centrifugal barrier and the short range of the nuclear force. However, when considering three nucleon scattering at 
higher energies, it appears natural to avoid a partial wave representation completely and work directly with vector 
variables. Only recently exact Faddeev calculations for three-body scattering in the intermediate energy regime be- 
came available. The formulation and numerical realization based on vector variables for the nonrelativistic Faddeev 
equations [13, [iMl as well as fully Poincare invariant ones [1^, [201 have been carried out for scalar interactions up to 
projectile energies of 2 GeV. 

Despite the technical sophistication with which the Faddeev equations in the continuum are solved today, the 
treatment of the singularity structure of the free three-nucleon propagator experienced only minor modifications from 
the original suggestion [ll| . e.g. in Ref. [l7[ the logarithmic singularities are integrated semi-analytically with splines 
in contrast to the earlier subtraction techniques. However, it would be most desirable to have a kernel without any 
logarithmic singularities. In Ref. [2l| a solution to this long lasting technical challenge is proposed and successfully 
carried out in the context of the partial wave decomposed Faddeev equations. It is the purpose of this paper to 
examine this suggestion in the context of a three-dimensional treatment of the Faddeev equation and introduce a 



kernel in which only simple poles in one variable occur. Those poles can than be integrated by standard subtraction 
techniques. 

In Section II we briefly revisit the form of the nonrelativistic Faddeev equation used in previous work [l7| to allow 
an easy comparison of the differences in our new approach. In Section III we describe in detail the simplification 
of the singularity structure of the free three-body propagator, and in Section IV we complete the calculation with 
the remaining angular integrations and connect to previous work in calculating the operators for elastic and breakup 
scattering. We conclude in Section V. 



II. THE FADDEEV EQUATION FOR THREE IDENTICAL BOSONS 



There are various presentations of three-body scattering in the Faddeev scheme y presented in the literature p, 
\12i . |23 | . We consider here the Faddeev equation for identical particles in the form 



T\cb)^tPl 



tPGoTl 



(2.1) 



The driving term of this integral equation consists of a two-body t-matrix t, the sum P of a cyclic and anticyclic 
permutation of three identical particles, and the initial state \4>) ~ jiys^qo), composed of a two-body bound state and 
the momentum cigenstate of the projectile particle. The kernel of Eq. (|2.ip contains the free three-body propagator, 
Go = {E — Hq + ie)~^, where E is the total energy in the center-of-momentum (cm.) frame. The operator T 
determines both the full breakup amplitude 



C/o - (1 + P)T 



and the amplitude for elastic scattering 



U = PGo 1 



PT 



(2.2) 



(2.3) 



For the explicit solution of Eq. ()2.ip the standard Jacobi momenta p, the relative momentum in the subsystem, and 
q, the relative momentum of the spectator to the subsystem are introduced. The momentum states are normalized 
according to (p'q'|pq) = S^{p' — p) S^{q' — q). Projecting Eq. (|2.ip onto these basis states leads to 



r(p,q;qo) = ro(p,q;qo) 

+ J d'p'd'p"d^q"tip, p'; e)(p'q|P|p"q 



-r(p",q";qo) 



(2.4) 



Here we abbreviate {pq\T\(pdqo) = r(p, q;qo) and the driving term as {pq\tP\ipd(lo) = Po(Piq;qo)- Under the 
integral we take advantage of the fact that the two-body t-matrix only depends on the relative momenta p and p' 
of the subsystem, and the thus the i-matrix is evaluated at the energy e = E — -^q^- The energy E" of the free 
three-body propagator is given by 



E"^-{p"'+'^q"') 
m 4 



The permutation operator is explicitly given as 

{p'q\P\p"c^')^S{p' + n,)S{p" 
with the 'shifted' momenta 



TTa) + 6{p' - TTl) S{p" + TTa) 



(2.5) 



(2.6) 



1 
TTl = -q 
^ 2^ 



T^2 



q+2q 



(2.7) 



Inserting Eqs. 



ai)- (P?7)) into Eq. ^^ leads to 

T{p, q, qo) = To(p, q, qo) + J d^p'd^p"d^q" t(p, p '; e) 

S{p' + TTl) ,5(p" - 772) + S{p' - TTl) 6{p " 



7r2j 



1 



T(p",q",qo) 



E + ie- E" 
The direct evaluation of the matrix elements of the permutation operator gives 



T'(p,q,qo) 



ro(p,q,qo) 
i(p,-7ri;£ 



(fq" 
1 



-T(7r2,q",qo) 



Defining a symmetrized t-matrix 



E + ie- E'' 

i(p, TTi; e)^7— ^ — r(-7r2, q", qo) 

E + le — E" 



tsip,7Ti;e) = t{p,7Ti;s) + f(p,-7ri;e) 



and realizing that for identical bosons T(7r2,q";qo) = r(— 7r2,q"; qo) one arrives at the expression 



7'(p,q,qo) 



To(p,q,qo) 

d^q"t,{p,TTi;e) 



E 



rn ^ -* 



r,ll2 



q q' 



-r(7r2,q",qo), 



3 

(2.8) 



(2.9) 



(2.10) 



(2.11) 



which is the starting point for the numerical calculations presented in Ref. [l7[. The free propagator in Eq. (|2.1ip 
clearly displays the difficulties inherent in numerically solving the three-body scattering problem in this form. The 
propagator depends on the magnitude of q" and through the scalar product q • q" on the angle between q" and a 
fixed axis given by q. The integration over cPq" leads to singularities with respect to that angle for each fixed value of 
q" and q. These singularities are integrable, but lead to logarithmic singularities in the variable q" . However, despite 
being integrable, they pose numerical challenges. For scattering calculations in a three-dimensional approach these 
challenges were met in Ref. [13, 13 for the nonrelativistic Faddeev equation and in Ref. [l^, [20| for Poincare invariant 
Faddeev equation. Nonetheless, three-body scattering calculations would be less challenging, if these singularities in 
the angles could be extricated from the ones in the momenta. A suggestion for accomplishing this task is made in the 
section III, following |2l| . 



III. A NEW LOOK AT THE SINGULARITY STRUCTURE OF THE FREE THREE-BODY 

PROPAGATOR 



As demonstrated in the previous section, using the vector variables in the evaluation of the matrix elements of 
the permutation operator of Eq. (|2.6p leads to the representation of the free three-body propagator of Eq. (|2.1ip . 
containing both integration variables q" and q • q". If we want to disentangle those two variables, we should start by 
separating the angular piece of the delta functions from the piece containing magnitudes of momenta and write 



This leads to 



Sip' + .,)Sip"-..) = ^^^^^^^^(p' + -r)^(p"--.) 
6ip'-.,)S{p" + n,) = '^'''j/'^'^''";j'h ip'-i.,)6ip" + i.,) 



p.. 



p" 



r(p,q,qo) 



Mp, q, qo) + f dp'p'^dp' dp"p"^dp" d^q"t{p, p', e 

6{p' + ivi) S{p" - ^2) + ^(p' - ^1) 6{p" + iv2 

S{p' - TTl) Sip" - TT2) 1 



E + ie- E' 



T(p",q",qo) 



= ro(p, q, qo) + / d'^q" dp' dp" 6{p' - ^1) 5{p" - ^2) 



E 



le 



±(-n"2 



(p"2 + |g"2) 



-(t,(p,p'7^i;e)r(p"7t2,q",qo) 



(3.1) 



(3.2) 



We arrived at the last expression by integrating over the angles and taking into account the symmetrized two-body 
t-matrix from Eq. (|2.10p . The magnitudes of the two vectors tti and 7r2 are given by 



iTTil = ^-q^ + q"^ + qq"x" 

k2| - ^q^ + p"^ + qq"x" (3.3) 

where cc" = q • q" . The above relations can be used to rewrite the two delta functions of Eq. (|3.2p in a form better 
suited for our further considerations. We rewrite 

Hp -^i)^% 5{x" - xo) 6(1 - |xo|) (3.4) 



where 



and 



'- I a 1 2 //2 1 ^ I 112 1 //2 2 \ fo c;^ 



8{p' - 712) = <5 ( P" - \Ip'' + p' - p"' elp"- Jp'^ + p^ - p"^ . (3.6) 



Inserting these expressions into Eq. ()3.2p gives 

r(p,q,qo) = To(p,q,qo) 

+ Jd'q" dp' dp"% 5{x" - xo) 6(1 - \xo\) 

5 L" ^p.2+ 3^2 „ 3^.2'! e [p'^ + Jg2 _ 3^"2 

E + ^e - l-}p"^ + y"^) ^^(P'^'^^ = ^) r(p"^2,q",qo). (3.7) 

Before we continue, we remember that the underlying two-body force supports one bound state with energy E^. 
This means, that is(p,p';z) has a pole at z = Ed- Because the transition operator T of Eq. ()3.7p is needed for all 
values of q, one will encounter this pole of t^ . Extracting the residue explicitly by defining 

, / / X is{p,p';z) 

t,(p,p;z) = — (3.8) 

Z ~ -Cyd 



and similarly for T and Tq one can rewrite Eq. p.7p as 
f(p,q,qo) = fo(p,q,qo) 



+ / d^q" dp' dp"% 5{x" - xo) 6(1 - l^ol) 

5 L" yp'2 + Jg2_3 A ^ (^^„ ^ 3^, _ 3^„, 

1 7 1 /^ , f(y^7r2,q",qo) ,„_. 

— -— W „2 I 3 „2\ ^s[p,PT^i\e) — — = 3-^;^. (3.9) 

Now we carry out the integration in p" and arrive at 

f(p,q,qo) = fo(p,q,qo) 

+ ^ y dq" dq" dp' S{x" - xo) 6(1 - \xo\) 6 ("p'^ + 3^2 _ 3^, 

A" r, ,. ^ f(p'^7r2,q",qo) 

-i,(p,p7ri;£) — — 3—-, (3.10) 



E + ze-Mp" + h')'' ' E + ie-Ed~i-q"-^'' 



where p" = \ p'"^ + f 9^ ~ 1 9"^- The two theta- functions in Eq. (|3.10p restrict the integration in the magnitudes of 

q" and p' into an area whose size depends on the magnitude of the spectator momentum q. An example of this area 
is indicated in Fig. [1] for the choice of g = 40 MeV/c. 

Next we consider the product of propagators in Eq. (|3.10p and separate this product as in Ref. [2l| 



1 



1 



E + ie- i(p'2 + |q2) E + ie-Ed- ^X'' 
1 1 



1 ('^'2 I 3„2\ 



:iP' 



ie- Ed 



3 „2 

4ni ^ 



-E, 



i,Q"' + iip" + h') 



The new denominator function, 



Giq,q",p'\ 



-E,-i^q"^ + l,ip'-2 + lq^) 



can not become singular inside the integration domain p' — q" . Using the relation from Eq. 

1 1 



we see that 



Giq,q",p'[ 



E, + ^v"'' \E4 + i:P' 



Lti"2 



>0. 



(3.11) 



(3.12) 



(3.13) 



With this separation of propagators the three-body transition amplitude from Eq. (|3.10p consists now of two integrals 
and can be written as 



r(p,q,qo) = To(p,q,qo) 

I r , , , 1 /•9/2+P' 

^ ~qJo ^^ E + te-^^{p'^+-lq^)J\,;2-p'\ 



dq"q"Giq,q",p') 



dq" 6{x" -xo)ts{p,p'7Ti;e) T{p" 7r2, q" ,qo) 
9 /-oo 1 rif^+i" 



dq"q" 

q Jo E + ie- Ed- 



_3_„//2 ,, , 
4my "'l9/2-9"| 



dp'p'G{q,q",p') 



f dc^'Six" ~xo)isip,p'TVi;e) f (p"7r2,q",qo) 



(3.14) 



In the second part of the kernel we changed the sequence of integrations over p' and q" . 

This new form of the kernel now exhibits only simple poles in thep' and q" integration, and the angle integration does 
not contain any singularity any more. As a remark, the simultaneous occurence of poles in the angle and momentum 
integration when solving Eq. (|2.1ip leads to the logarithmic singularities, which are not present in Eq. (|3.14p . The 
poles of the free propagator occur in the variable p' and the deuteron pole in the variable q". 



IV. THE ANGLE INTEGRATION 



Since we ignore spin and isospin dependencies, the matrix element T(p, q, qo) is a scalar function of the variables 
p and q for a given projectile momentum qg. As was shown in Ref. [l7[, one needs 5 variables to uniquely specify 
the geometry of those three vectors. For the clarity of presentation we repeat some of the arguments here. Having 
in mind that with three vectors one can span two planes, i.e. the p-qg-plane and the q-qg-plane, a natural choice of 
independent variables is 



P 



IPI: Q 



P ■ qo, Xg^q- qo, xl" = (qo x q) • (qo x p) 



(4.1) 



The last variable, x^^, is the angle between the two normal vectors of the p-qo-plane and the q-qo-plane. It should 
further be pointed out, that the angle between the vectors p and q, ijpq = p • q is not an independent variable, but 
can by related to the ones given above as 



For the special case where qo is parallel to the z-axis ((/o-system) one can write 

Vpq = XpXq + y 1 - X^Jl - X\ COS^Ppq, 



(4.2) 



(4.3) 



where the angle fpq is the difference of the azimuthal angles of p and q related to the specified z-axis. 

The (5-function in the angle integration of Eq. (|3.14p . 5{x" — Xq) = 5{i\ ■ q" — Xg) suggests to choose for the q" 
integration the z-axis parallel to the vector q (g-system). The angle dependence of the two-body t-matrix is then 
explicitly given as 

4(p,p'7ri;-E(g)) = 4b,p',p-7ri;e), (4.4) 

where 

P ■ ^1 = - ^p'^ + y^'^" (4.5) 



and 



Vpq" = P • q = VpqX- 



" + ^1 - Vlq Vl - x"2 cos(^p - ^"). (4.6) 



The angle Lpp is the azimuthal angle of p in the g-system. As mentioned before, x" = q • q" and q" = |q"|. 
The angle dependence of the three-body transition amplitude is more intricate and can be written as 

f(p"7r2,q",qo)=r(p",x,,,x* „ ,y,,,„ ,q";qo). (4.7) 

Explicitly these angles arc given as 

92;^ + 2? Uqoq" 

7r2 • qo = Xtt^ - 



f + W'^ + qq"x" 



q"qo 



Vqoq" = XqX" + ^ 1 - X^ \J 1 - x""^ COs{Lpq„ - (/?") 



qx ~^q 

,0 ^ ^2-g"-(^2-go)(g"-go) ^ ^q^+y^+qq-'^- '^"^^^°^" .^g. 

"^^^-" VI - (-2 • ^o)^ VI - iq" ■ qor V^^^../^^<7' 

Like ipp in Eq. ()4.6p the angle (pg^ is the azimuthal angle of cfo in the g-system. It was shown in Ref. [l7| that because 
of the if" integration, only the knowledge of cos{ipqg — (pp) is required. This difference can be explicitly represented as 

cos(v.,„ - ^p) = /"■P-(g-^°)(g-^)^ = -P - ^^^"^ . (4.9) 



v/1 - (9 • 9o) Vl - (P ■ 9)' \A^y^ 



l-x2 



Since only difference of the angles enters, one can choose ipq^ arbitrarily, e.g. ipq^ = 0. Furthermore, cosiy9p and simpp 
required in Eq. (|4.6[) arc then also uniquely given [17| . 

With these preparations we are ready to carry out the angle integration J dq" = J_i dx" L" dip" of Eq. p.l4p 
explicitly. The x" integration is fixed by the (5-function in terms of xq = xo{q,p' ,q") from Eq. (|3.5|) . leaving only an 
integration over ip" . Explicitly, the variables of Eqs. (|4.6p and (|4.8p need only be evaluated at a fixed x" = xa{q,p', q"). 
Thus, the explicit representation for the transition amplitude T reads 

f{p,Xp,xl°q,Xq,q;qo) = fo{p,Xp,xl°q,Xq,q;qo) (4.10) 

+ - dp'p'——, 1 , ,„ 3 „, / dq"q" G{q,q",p') 



I 



dip ts p,p ,— f — , 

° V ^y^+q"'^ + qq"xn 



T\P , I ,xl\y^^^„(xo),yqoq"[xo),q ; go 



2 r 
q Jo 



1 rq/^+q" 

^^"^"tt^ 7^ ^^ / rfpV G(g, g",p') 



2ir 



dip ts\p,p,- 



T I r," g^g + 2l"yQol"i^o) qo (^ \ ,, „^„ \ „".„ 

/q^ + y^ + qqoxo 

where p" = wp'^ + f?^ — |(Z"^ is fixed. 

The only remaining detail is to provide an explicit expression for the Born term, 

r(p,q,qo) = (pq|tP|</Pdqo), (4.11) 

where (pd stands for the deuteron bound state. Projecting on Jacobi momenta and evaluating the permutation operator 
leads to 

r(p,q,qo) = '^d(q+7:qo )is ( P, 7:q + qo, e ) ■ (4.12) 



2^"y V' 2 
Using the invariant variables of Eq. (14. ip one arrives at 

fo{p,Xp,xll,Xq,q]qo) = (pd { \l q^ + 7Q0 + IQa^q ) (4.13) 




The transition amplitude of Eq. ()4.10|) together with the Born term given above has as far as the angle integration 
is concerned a structure very similar to the one given in Eq. (2.19) of Ref. [l3|. Even in the case where logarithmic 
singularities have to be integrated, they are independent of the angle Lp" . Thus, in the numerical realization demon- 



strated in Ref. [17[ the (^"-integration of the kernel is carried out for each fixed value x" and q" on their respective 
grids. Then the singularity structure depending only on the x" and q" variables is explicitly dealt with. In Eq. (|4.10p 
the (/^''-integration needs to be carried out for each fixed value q" and p' given on their respective grids. These grids 
also fix xo{q,p',q"). Since the functional dependence on ip" is the same in both integrals over ip" as is the area of 
integration in the p'-q"-plaTie, the integral over ip" in the kernel needs to be evaluated only once if both integrals 
of Eq. (|4.10p are calculated on the same p'-q"-giid. Under these conditions the numerical effort as far as the angle 
integration is concerned is similar to the one of Ref. [131 • 

Once the transition operator T(p,q, qo) is explicitly calculated as function of the 5 independent variables, the 
amplitude for elastic scattering is obtained by calculating the matrix elements of the operator U given in Eq. (|2.2p as 

{q(pd\U\qQipd) ^ 2(/3d f-q + qoj f £■- — (g^ + q-qo + go)) "y^d (q+ 2^0 



2 / d^q"p,d -q + q" ^ % , ^. (4.14) 



The amplitude for the full breakup process is given according to Eq. (j2.3p by 

(pq|f|qo¥5d, 
{pq\Uo\qo(pd) = 



E^^^q^-Ed 

\P + |q, -P - iq|r|qo(^d) 1-\V - |q, +P - \(i\f\q_aPd 



E-i,{~V-W-Ed E-^{+p-\qY-Ed 



(4.15) 



Both operators are explicitly given in Ref. [17| using the independent variables of Eq. (|4.ip . and can be directly applied 
with the expression of Eq. (|4.10p . 



V. SUMMARY AND CONCLUSIONS 

In Ref. [iTJ l the formulation of the nonrelativistic Faddeev equation for three identical bosons as function of vector 
variables was introduced and successfully solved for laboratory projectile energies up to the GeV regime. The key 
point allowing the calculation at those higher energies is to neglect the partial-wave decomposition generally used 
at lower energies and to work directly with momentum vectors, thus including all partial waves automatically. In 
the formulation of the Faddeev integral equation in the continuum which is most widely used, the singularities of 
the free three-body propagator occur simultaneously in an angle and momentum integration, leading to the socalled 
logarithmic singularities. They require special care in numerical applications. Although sophisticated algorithms have 
been developed to integrate those singularities along the real axis, it still is desirable to have a formulation of the 
Faddeev kernel, in which the singularity structure is sirnpler. 

Starting from the formulation given in Refs. [I4I and l2l| we propose a new formulation of the Faddeev kernel which 
does not contain this technical obstacle of logarithmic singularities. Instead, the singularities of the free threc-nucleon 
propagator appear as poles in a single variable. Those kind of poles can be integrated with standard techniques as 
e.g. used in integrating the two-body Lippmann-Schwinger equation. The singularity given by the deuteron pole 
is a simple pole (as before), but is now cleanly separated in a separate integration. The integration over the angle 
variable not affected by the pole structure remains very similar. This simplification in handling the singularities of 
the three-body continuum in a similar fashion as the two-body continuum should ease applications of the Faddeev 
integral equations in different areas of physics. 
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FIG. 1: (Color online) The domain for the integration over the momenta p' and q" as function of the external momentum q. 
Here q — 40 MeV/c is chosen. The area of integration is the rectangle enclosed by the three lines. 



